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ABSTRACT 



The very first stars in the Universe can be very massive, up to IO^Mq. If born 
in large numbers, such massive stars can have strong impact on the subsequent star 
formation, producing strong ionising radiation and contaminating the primordial gas 
with heavy elements. They would leave behind massive black holes that could act 
as seeds for growing supermassive black holes of active galactic nuclei. Given the 
anticipated fast rotation, such stars would end their live as supermassive collapsars and 
drive powerful magnetically-dominated jets. In this letter we investigate the possibility 
of observing the bursts of high-energy emission similar to the Long Gamma Ray Bursts 
associated with normal collapsars. We show that during the collapse of supercollapsars, 
the Blandford-Znajek mechanism can produce jets as powerful as fewx lO^^erg/s and 
release up to lO^^erg of the black hole rotational energy. Due to the higher intrinsic 
time scale and higher redshift the initial bright phase of the burst can last for about 
10^ seconds, whereas the central engine would remain active for about one day. Due to 
the high redshift the burst spectrum is expected to be soft, with the spectral energy 
distribution peaking at around 20-30keV. The peak total flux density is relatively low, 
10~^ergcm^^s~^, but not prohibitive. If one supercoUapsar is produced per every mini- 
halo of dark matter arising from the 3-tT cosmological fluctuations then the whole sky 
frequency of such bursts could reach several tens per year. 

Key words: cosmology: early Universe - black hole physics - gamma-rays: bursts - 
X-rays: bursts - radiation mechanisms: non-thermal - relativity - magnetic fields 



1 INTRODUCTION and the eventual accumulation of mass proceeds via ac- 
cretion of the surrounding gas. Very high accretion rate, 
Accordmg to the modern hierarchical clustermg theories oi . ,„_3,, _i i. .. .i n i . r 

, ,^ , , „ , M>Mc~4xlO M0yr , may limit the final mass to few 

t , ■ galaxy formation the first stars are born withm collapsed i i i i j_i ^ ^ n i • ^ 

^f, _,, . , hundreds solar masses as the protostellar luminosity reaches 

Ci ■ haloes of dark matter of 10 at z ^ 20 The pri- Eddington limit JOmukai fc Pallall2003l l. For lower ac- 

mordial gas falls into the potentiaf weff of these hafoes and , . , . , , . , n . i 

, cretion rate, the accretion may proceed even alter the onset 

fragments into clamps of ~ 10 Mr via gravitational msta- r i i • • ^ n j • 

, r— — -fi -1 „ , , „ or nuclear burning m the stellar core and result m the fi- 

bility (IBromm et al.l 120021 ) . Because this gas is metal- free , , , i ^3 » <■ ivt • i j_ i • r i • i 

* ; — — " . . . _ nal mass Al ~ 10 AIq. Numerical studies ot cosmological 

its cooling IS rather slow and further fragmentation into ^. , . j_ , -i-j, j_ j_i ^ i • • • i 

„ , ^ , . , , / r I rrV ^ \ — rn u^^nl gTavitatiouaf instability suggest that, although m principle 

smaller clamps seems to be avoided let. I lurk et al. M2009I : .. , , i • i r 10-2^^ -1 • 

r- , x T the accretion rate can be as high as lewxlO Mqyv , m 

Ibtacv et al. 2009). Instead, the clumps contract m a quasi- 1.. .1 ... 1 . ■ , .. n 1 

~; — , . ' , . , , r. reality the rotational support against gra vity often becom e 

static tashion as a whole, suggesting that the first stars can . ^ j_ 1 j 1 1 »V J7T 1 — nUr,r,JN 

, ^ rr^n , . . . , importaut and reduces the rate below Me UGao et al.ll2007l ). 

be very massive indeed, M > lOOM^. ihe actual initial [Trr; — ; ; — n JT^TJT^ j_ j- 1 1 i- r j_- n 

/ . ' , lUhkubo et al.l II2009I I studied the evolution ot accreting Pop- 
mass tunction (IMl ) or first metal- tree stars Population 1 . • ttt . <■ ,1 ■ 1 , 

^ ' r ulation 111 stars trom tfie pre-mam sequence evolution to 

111 stars), however, is not known yet as too many factors ,, , , .1 . n 1 1 

' the core-collapse and conhrmed that the hnal mass can be 
come into play, making the problem intractable analytically, , i^3j^-ir ■ n ^ ^ ^ 

/ " T s-s large as 10 Mq. Very massive farst stars are also pre- 

and rather challenging numerically. In particularly, the mi- . 1 . .1 . . 1 . 11 -i -i , . / 

„ _o.. dieted m theories inv olving dark matter annihilation (e.g. 

tial mass ot protostars can be very small, down to 10 Mq, iNataraian et allbood ) 

Population III stars with masses I4OM0 < M < 

* serguei@maths.leeds.ac.uk (SSK) 26OM0 most hkely end their life as pair-insta bility super- 

t E-Mail: bmv@maths.leeds.ac.uk (MVB) novae which leave no compact remnant behind (|Fryer et al.l 
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I2001I ). If such stars were the main outcome of initial 
star formation they would overproduce heavy elements 
in the early Universe, in conflict with the observations 
of extremely metal pure stars in galactic bulges and 
the obser ved abundances of intergalactic and intercluster 
medium (lUmeda fc Nomotd I2OO2I : iHeeer fc WooslevI I2OO2I : 
IChieffi fc Limongill2002l ). 

More massive stars, which will be referred to as Very 
Massive Stars (VMSs) are expect ed to collapse int o black 
holes with very little mass loss l|Frver et al.l I2OO1I ). They 
would leave behind massive black holes (MBHs), which 
could play the role of seeds for the supermassive black holes 
(SMBHs) of Active Galactic Nuclei (AGNs). Assuming that 
MBHs are formed at t he rate of one per min i-halo developed 
from a 3(T-fluctuation. lMadau fc Reed l|200ll ) estimated their 
density to be around 5000 per galaxy like Milky Way and 
their total mass comparable to the total mass of SMBHs. 
This suggests that SMBHs could form via mergers of MBHs, 
the idea that has being actively developed in recent years. 
Even more massive VMSs, with M ~ 10® Mq, could be 
formed in more massive dark matter ha loes, with total 
mass M ~ 10^ Mq, collap sed at z ~ 10 (|Bromm fc Loebl 
l2003l : lB"egelman et al.ll2006l ). Although much more rare, such 
events can provide an alternative way of producing SMBHs. 

From the observational perspective it is difficult to dis- 
tinguish between a VMS and a cluster of less massive Pop- 
ulation HI stars. This suggests to investigate the potential 
observational consequences of VMS collapse, which could 
be quite spectacular because of the very high mass involved. 
Given their expected fast rotation, it seems li kely that su- 
ercoU apsars (classified as type-Ill coUapsars in lHeger et al.l 
20031 )) develop accretion disks, drive relativistic jets, and 
produce bursts of high energy emis sion in the fashion sim- 
ilar to their le ss massive relatives ("MacFadv en fc WooslevI 
I1999I : [Barkov fc Komissarov 2008). If detected, they would 
become the most distant strong sources of light and pro- 
vide us with a new way of probing the physical conditions 
in the epoch close to the "Dark ages". Even non-detection 
could be useful, allowing to set constraints on models of star 
formation in the early Universe and the origin of SMBHs. 
There has been a number of papers looking into cosmolog- 
ical evolution of GRBs, including those from the Popula- 
tion III stars fe.g.lBromm fc LoebllioO^ : iLe fc Dermej|2007l : 
iNaoz fc Bromberell2007l ). They assumed that the Population 
HI GRBs are similar to those from the lower redshift Pop- 
ulation II stars whose mass is significantly below IOOM0. 
However, because of their very high mass and redshift, the 
GRB-like bursts of supercoUapsars can be rather special. 

There are two crucial differences between a normal col- 
lapsar and a supercoUapsar. One is that instead of a proto- 
neutron star of solar mass the supercoUapsars develop proto- 
black holes of tens of solar masses, within which the neutri- 
nos from e l ectron capture are trapped (Frvcr et al. 2001; 
ISuwa et al.l bOOTh . The other is that the accretion disks 
of supercoUapsars are far too large and cool for the neu- 
trino annihilation mechanism. This has already been seen 
in the numerical simulations of supercoUapsar with mass 
M = SOOMq (Hrycr et al. 2 00il). Ut i lising the study of 
hyper-accreting disks by iBeloborodovl (|2008l ) we find that 
at best the rate of heating due to this mechanism is 



where M is the accretion rate and Mh is the black hole mass. 
(Here and in other numerical estimates below we use the fol- 
lowing notation: is the mass accretion rate measured in 
the units of IO^Mqs"^ and Mk is the mass measured in 
the un its of 10*° M©.) Such low values have lead lFrver et al.l 
l|200ll ) to conclude that the magnetic mechanism is the only 
candidate for producing GRB jets from supercoUapsars. In 
the following we analyse one particular version of the mech- 
anism, namely the one where the jets are powered by the ro- 
tational energy of the black hole v i a the Blandford-Znajek 
proce ss (|Blandford fc Znaiekl 1 19771 : iBarkov fc Komissarov I 
I2OO8D . 



2 BLANDFORD-ZNAJEK JETS FROM 
SUPERCOLLAPSARS 

VMSs are expected to rotate rapidly, close to the break-up 
speed and hence produce rapidly rotating MBHs. Moreover, 
in the absence of strong magnetic field in the prestine pri- 
mordial gas, VMSs will be weakly magnet ized and as the re- 
sult c ould develop rapidly rotating cores (jWooslev fc Hegerl 
I2OO6I ). This suggests that the spin parameter of MBHs can 
be very high, a ~ 1, yeilding the enormous rotational energy 
_Erot — 5 X 10^®Mh,3erg. In order to estimate the Blandford- 
Znajek luminosity we need to know the strength of magnetic 
field accumulated by the hole. The usual approach is to re- 
late it to the gas pressure in the disk and for this we need 
to know the parameters of the accretion disk itself. Accu- 
rate determination of these parameters and their time evo- 
lution requires to know the structure of VMS prior to the 
collapse and the physics of accretion disk. Unfortunately, 
this information is lacking at the moment. In particular, al- 
though the structure of supermassive stars. Ma ^ 10^ M© 
is very well described by a polytropic model with n = 3 
(|Zeldovich fc Novikovlll971 ), we are mor e interested in stars 
with < lO^MfTi. iFrver et all l|200ll ) studied the struc- 
ture and evolution of a 3OOM0 star. Prior to the collapse 



£ ~ 2 X lO"*** Mn''''M, 



-3/2 

h,3 



erg s 



(1) 



the star entered the red giant phase and expanded from 
i?s = 4 X lO^^cm to 1.5 X lO" cm. However, the initial ro- 
tation rate of this star was slow compared to about ~ 50% 
of the brake-up speed invoked the current single-star model 
of GRB progenitors. In this model, the progenitors remain 
chemically h omogeneous and compact all the way up t o 
the collapse (|Yoon fc Langeij|2005l : IWooslev fc Hegej|2006l ) . 
Given this lack of i nformat i on ab out the progenitor struc- 
ture we will follow iBethd l|l990l ) and assume p (x R 
{Rc <R <Rs) distribution of mass density in the star prior 
to collapse. In fact, this distribution agrees reasonably well 
with the numerical models of rapidly rotating low metallic- 
ity stars considered as likely progenitors of normal GRB (see 
Figure 2 in Kumar et a l. 2008). As to the stellar rotation, 
we will assume that it is uniform (f2=const) in the stellar 
envelope, with 50% of break-up speed at the stellar surface. 

Due to the slow neutrino cooling, the accretion disks 
of supercoUapsars are expected to be radiatively inefBcient, 
with possible exception only for the very inner region. This 
suggests to us e the ADAF (Advect ion Dominated Accretion 
Flow) model (|Naravan fc Yilll994l ) to describe these disks. 
Since the radiation pressure dominates, we can use the ratio 
of specific heats 7 = 4/3 which gives us 



© 0000 RAS, MNRAS 000, 000-000 



SupercoUapsars and their X-ray Bursts 3 



3a 



2 2 
— Vh, 

7 ' 



H ~ Rcs/vk, 



(2) 



where is th e accretion speed, Cs is the sound speed, 
Vk = \J GMjR is the Keplerian speed, H is the vertical 
disk scale, and a is the effective viscous stress parameter 
of the Q-disk model (jShakura fc Sunvaev I [l973 ) . The disk 
density and pressure can be estimated combining the above 
equations with the expression for the mass accretion rate, 
M ~ 2nRH pvin- Straightforward calculations yield 



(3) 



14 MjGMf'^ 
12-Ka ^575 



The poloidal magnetic field should scale with the thermo- 
dynamic pressure, so we write = SnP/P, where /3 is the 
magnetization parameter. Applying this equation at the ra- 
dius of the marginally bound orbit -R„b = fi{a)Rg, where 
fi{a) = 2-a + 2(1 - a)^^'^ and Rg = GMh/c^ is the BH's 
gravitational radius, we find 

B,„b ~ 3 X lo"/r'^''''/3y^A^o^^Gauss, (4) 

where /3i = /3/1O0 If the magnetic field is generated in 
the disk then it is likely to change polarity on the scale 
~ H. This may lead to significant variation in the strength 
and polarity of the magnetic field accumulated by the 
black hole and reduce the Blandford-Znajek luminosity (e.g. 
iBarkov fc Bau shcv''200^. We will assume that this effect is 
accounted for in the value of /3. 

The power of jet energized via the Blandford- 
Znajek mechanism can be estimated using the monopole 
solution for magnetosph eres of rotating black holes 
l|Blandford fc Znaiek|[l977l ). which gives 



3c V 47r / 



(5) 



where Qh = f2{i)(? /GMh is the angular velocity of 
BH, /2(a) = a/2(l -f \/l — a^), and is the magnetic 
fiux threading one of the BH's hemispheres. Inside the 
marginally bound orbit the disk plasma quickly dives into 
the BH and the mag netic fiux c an be roughly estimate as 
* = 27r_Rib5„,b ( Reynolds et al.ii2006. '). Combining this re- 
sult with Eas. (l4l5p . we find 



14 fl'^f. 



9 a/3 



Mc 



0.05 



a-i 



/3i 



Mc 



(6) 



where a_i = a/0.1 ( for 0.5 < a< 1 the combination ./f'^^/l 
depends weakly on a and is approximately 1/4). 

The mass accretion rate can be estimated f ollowi ng the 
procedure described in iBarkov fc Komissarovl l|2009l ). The 
total accretion time includes the travel time of the rarefac- 
tion wave send into the stellar envelope by the core col- 
lapse, the time of the envelope collapse, and the disk ac- 
cretion time, which gives the largest contribution. Account- 
ing only for the disk contribution, the accretion time scale 
for the stellar matter located in the progenitor at radius 
R is t <x f IM^, where I = Q.R^ and M(R) is the stel- 
lar mass enclosed within radius R. Then M = dM/dt ~ 
[dM / dR) / [dt / dR) , where for the Bethe's model we have 



^ The inner edge of the disk varies between the marginally bound 
and the marginally stable orbits, depending of its thickness. In our 
calculations this does not make much difference. 



dM/dR ~ Ms/{R\n{Rs/Rc)) and dt/dR ~ Qt/R, where Rc 
is the stellar core radius. Collecting the results we obtain 



M 



: 36- 



6ln{Rs/Rc) t t ' ^'^^ 

where t is measured in seconds and we used Rs/Rc = 1000 
Here, we assumed that the whole of the disk is accreted by 
BH, following the original ADAF model. However, it has 
been argued that this model has to be modified via includ- 
ing disk wind (Advection Dominated Inflow Outflow Solu- 
tion (ADIOSl. fBlandford fc Begelman|[l999t ). which implies 
a mass loss from the disk and a smaller accretion rate com- 
pared to Eq.(l7]). While the arguments for disk wind are very 
convincing, the actual value of mass loss is not well con- 
strained and can be rather low. Given Eq.Q the power of 
BZ-jet is 



3.2 X 10""e„ Ms,3 -1 
-ergs , 



a_i/3i 



t2 



(8) 



where t2 = t/100 and < 1 is the fraction of the disk mass 
reaching the BH. Given the jet propagation speed inside the 
star, Vj 0.2c, deduced from axisym metric numerical simu- 
lations l|Barkov fc Komissarov |[2008l ). the jet breakout time 
is expected to be around of few hundred seconds and, thus, 
the numerical factor in Eq.([8ll gives us the optimistic jet 
power at the time when it becomes observable. In fact, the 
initial influx of mass through the polar column is very large 
and acti vation of the Blandford-Zna jek mechanism can be 
delayed (jKomissarov fc Barkovll2009l ). The very latest time 
for the activation is given by the free-fall time of the whole 
star. 



1/2, 



(9) 



as by this time the polar column becomes completely empty. 

The total duration of the jet production phase has to 
be similar to the disk lifetime. If VMS is rotating at half of 
the break-up speed then the initial outer edge of the disk is 
at Rd ~ Rs/4:- Ignoring the edge expansion due to accumu- 
lation of angular momentum, the disk life time is given by 
its "viscous" time scale 



2Rs 



3«in(i?s) 



: 5000aZlRt%M, 



-1/2, 



(10) 



where Rs,i2 is the stellar radius measured in lO^^cm. By 
this time the BZ power will be significantly reduced but 
could still play a role in shaping the li ght curve of afterglow 
emission (|Barkov fc Komissarovll2009l '). 

Using the mass accretion rate given by Eq.lO we can 
check if the neutrino cooling needs to be included in the 
model. Under the conditions of the supercoUapsar's disk its 
cooling is dominated by pairs . Using the well know n equa- 
tion for this cooling rate (e.g. lYakovlev et al. Il200ll ) we can 
compare the cooling time with the accretion time at a given 
disk radius. The result is 



^ ~ 0.3a:?/*(7^/i^,)-l3/8M5/4M7|/^ 



(11) 



Thus, except for the very inner part of the disk, the neutrino 
cooling is indeed inefficient. 



2 For a 300Mq star at t ~ 8s this gives M ~ 1.3Mi7)S~l , whic h 
agrees reasonably well with numbers given in I Fryer et al.l 1I2OOII) . 
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The high BZ power given by Eq.((8]) suggests that the 
GRB-like burst emission from such jets could be seen even 
from 2; ~ 20 and in the next section we discuss the properties 
of such bursts in more details. 



3 OBSERVATIONAL SIGNATURES 

Assuming that the radiation mechanism of the supercoUap- 
sar jets is similar to that of normal GRB jets, we expect the 
peak in the spectral energy distribution of the prompt emis- 
sion in the source frame to be around O.SMeV. However, the 
cosmological redshift effect reduces the peak down to 



25keV 



+ 

V 20 / 



(12) 



which is still inside the energy window of Swift's BAT. For 
the same reason the observed total duration of the burst 
increases up to 

l + z\ 



tb 



20 



— 1 7->. 

a-iR 



1 r>3/2 



1/2 



day. 



(13) 



The characteristic source frame time scale for the decay of 
BZ luminosity in the model presented above is given by the 
time since the onset of the collapse. Thus, the initial time 
scale for the burst decay will be of order of the jet break out 
time, fewxlO^ seconds, or a bit longer if the activation of 
the BZ-mechanism is significantly delayed. In the observers 
frame this translates into fewx 10'^ — 10* seconds. Thus, these 
bursts would appear not only unusually soft but also uimsu- 
ally long-lasting. 

The total fiux density of the burst emission received 
on Earth and the isotropic luminosity are related via F = 
L/( 47rr^,), where T l is the luminosity distance to the source 
fe.g. |Peeblesll993l ). However, the emission from GRB jets is 
highly anisotropic due to the relativistic beaming. Moreover, 
not all of the Blandford-Znajek power is converted into the 
radiation within the energy window of the receiver. This 
leads to 



F : 



(14) 



where ^ ^ 1 is the solid angle of the radiation beam and 
and Ec < 1 is the conversion efficiency. In flat Universe 



Ho 



(i + z) + zf + Qa) 



-1/2 



dz. 



(15) 



due to internal shocks in baryon dominated flow, as this is 
proposed in the curr ently most popular mo del of prompt 
gamma-ray emission eszaros fc Rees I ll99i ). then the su- 
percoUapsar burst produced by magnetically-dominated BZ 
jet may well be smooth and featureless. However, there are 
models of normal GRBs that attribut e the observed variabil- 
ity to unsteady magnetic dissipation ( Lvutikov fc BlandfordI 
l2003l: lGiannios et al ."2009': Kuma r fc Narava n 2009). If they 
are correct then the supercoUapsar bursts will also show fine 
substructure. 

In order to esti mate the observ e d rat e of such burst 
we assume, following iMadau fc Reej (|200ll ). that the dark 
matter mini-haloes that host supercoUapsars arise from 3-a 
fluctuations that constitute only ~ 0.3% of the dark matter 
matter of the Universe and that only one supercoUapsar per 
mini-halo is produced. The total mass per Mpc^ at « = 20 
is Mmpc — 1.5 X 10^^ Mq. The number density of 3-a mini- 
haloes is then 

,nd,„MMpo 



: 0.003- 



lO'^Mpc"^ 



106 M© 

Let us assume, for the sake of simplicity, that all supercol- 
lapsars go off simultaneously at cosmological time te corre- 
sponding to 2 = 20 (a moderate spread around this redshift 
will not significantly change the result). In fiat Universe the 
observed time separation between events occurring simulta- 
neously at ro and ro + dro, where ro is the comoving radial 
coordinate, is dto = cdro. The corresponding physical vol- 
ume within one steradian of the BAT's field of view, is 

dV = a^{te)rldro, 

where a{tc) = (1 + z)~^ is the scaling factor of the Universe 
at t = te (in the calculations we fix the scaling factor via 
the condition a{to) = 1). ro and te are related via ro = 
rL(l + z)~^ . Putting all this together we find the rate to be 



fc^A 



cw„hrL 
(1 + 2)5 



iw) 



-1 -1 
yr sr 



Recent high-resolution simulations of cosmological star for- 
mation indicate the possibility of further fragmentation of 
gas clumps in minihaloes, resulting in formati on of binary 
or even multiple protost ars in some realizations ijTurk et al. I 
I2OO9I: IStacv et al.l l2009''l. Thus, the theoretical rate of VMS 
formation can be significantly smaller compared to the one 
used in our calculations, making the supercoUapsar bursts 
rare events. This may explain why such bursts have not been 
seen so far. 



For 2 = 20 and the densi ty parameters SIa = 0.72, Qm = 
0.28 l|Komatsu et al.ll2009l ') this gives us 

f ~ 2 X 10~^ ^"'~'^T M.,3t3~'ergcm-"s-\ (16) 

where 6^,-1 = tc/O.l, As = A/W'^, and t-j = t/10^. One 
can see that for the first 10*s this is above the sensitivity of 
BAT, 10~*ergcm~^s~^, and thus such a burst could trigger 
BAT. Having said this, we keep in mind that there is a 
great deal of uncertainty with respect to the values of various 
parameters appearing in Eq. (|16|l . 

The time dependence in Eq. lfT6)l gives the evolution of 
mean bolometric fiux. It is not clear if the supercoUapsar 
bursts will also exhibit the fine substructure characteris- 
tic of normal GRBs. If the variability of normal GRBs is 



4 CONCLUSIONS 

In spite of the significant progress in the astrophysics of 
Gamma Ray Bursts, both observational and theoretical, it 
may still take quite a while before we fully understand both 
the physics of the bursts and the nature of their progenitors. 
At the moment there are several competing theories and 
too many unknowns. Similarly, we know very little about 
the star formation in the early Universe. For this reason, 
the analysis presented above is rather speculative and the 
numbers it yields are not very reliable. Further efforts are 
required to develop a proper theory of supercoUapsars and 
to make firm conclusions on their observational impact. On 
the other hand, our estimates suggest that if we are on the 
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right track then the X-ray bursts of supercoUapsars may 
be detectable already with Swift. The expected very long 
duration of bursts and their relatively low brightness imply 
that a dedicated search program using the image trigger may 
be required. Such search would be useful even in the case 
of non-detection as this would put important constraints on 
models of star formation in the early Universe, models of 
the GRB progenitors, and the origin of SMBHs. 
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